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Abstract. It was shown recently [R. Behrend, P. Di Francesco and P. Zinn- Justin, On 
the weighted enumeration of alternating sign matrices and descending plane partitions, 
J. Combin. Theory Ser. A 119 (2012), 331-363] that, for any n, there is equality between 
the distributions of certain triplets of statistics on n x n alternating sign matrices (ASMs) 
and descending plane partitions (DPPs) with each part at most n. The statistics for an 
ASM A are the number of generalized inversions in A, the number of — I's in A and the 
number of O's to the left of the 1 in the first row of A, and the respective statistics for a 
DPP D are the number of nonspecial parts in D, the number of special parts in D and the 
number of n's in D. Here, the result is generalized to include a fourth statistic for each 
type of object, where this is the number of O's to the right of the 1 in the last row of an 
ASM, and the number of {n — l)'s plus the number of rows of length n — 1 in a DPP. This 
generalization is proved using the known equality of the three-statistic generating functions, 
together with relations which express each four-statistic generating function in terms of its 
three-statistic counterpart. These relations are obtained by applying the Desnanot-Jacobi 
identity to determinantal expressions for the generating functions, where the determinants 
arise from standard methods involving the six-vertex model with domain-wall boundary 
conditions for ASMs, and nonintersecting lattice paths for DPPs. 



1. Introduction 

Alternating sign matrices (ASMs) and descending plane partitions (DPPs) are combi- 
natorial objects which share some basic enumerative properties. (For an introduction to 
ASMs, DPPs and related subjects see, for example, the reviews of Bressoud [ ], Bressoud 
and Propp [ ], Propp [50], Robbins [ . or Zeilberger [G3]. For a more up-to-date, but briefer, 
introduction see [■), Sec. 1].) 

ASMs and DPPs initially appeared in different contexts, but at approximately the same 
time, and it was shortly thereafter that the enumerative connection between the two types 
of object was first observed. In particular, a conjecture was formulated by Mills, Robbins 
and Rumsey Conj. 3], stating that the distributions of certain triplets of statistics on 
n X n ASMs and on DPPs with each part at most n are equal for any positive integer n. 
Interestingly, at this time, an aspect of the unproved conjecture was helpful to Mills, Robbins 
and Rumsey in proving a different conjecture, namely the Macdonald conjecture for the 
generating function of cyclically symmetric plane partitions [ ' ']. 
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Some special cases of the conjecture relating ASMs and DPPs were subsequently confirmed 
(see [5, Sec. 1.3] for a detailed overview of these), and the full conjecture was proved recently 
by the present authors [5, Thm. 1]. Among the previously-proved special cases, probably 
the best-known is simply that the total numbers of n x n ASMs and DPPs with each part 
at most n are equal for any n. This follows from formulae of Andrews [ ] for DPPs, and of 
Zeilberger [60] or Kuperberg [3(S] for ASMs, which give the numbers of each type of object 



The primary aim of this paper is to generalize the result of [5, Thm. 1], involving equally- 
distributed triplets of statistics on ASMs and DPPs, to a result involving equally-distributed 
quadruplets of statistics. The three original statistics for an ASM A are the number of 
generalized inversions in A, the number of — I's in A and the number of O's to the left of 
the 1 in the first row of A, and the respective statistics for a DPP D with each part at 
most n are the number of nonspecial parts in D, the number of special parts in D and the 
number of n's in D. The additional statistics considered here are the number of O's to the 
right of the 1 in the last row of an ASM, and the number of (n — l)'s plus the number of 
rows of length n — 1 in an DPP with each part at most n. In each case, the first and second 
statistics depend on the bulk structure of the object, while the third and fourth statistics 
depend only on the structure at or near its boundary. 

The joint distribution of the fourth ASM statistic with some or all of the other three ASM 
statistics has previously been considered. See Section 7.3 for an overview of such studies 
in the literature. However, it seems that the fourth DPP statistic has not appeared before 
in the literature, and that no other previously-studied DPP statistic has the same joint 
distribution with the other three DPP statistics as that of this statistic. 

In studies of ASMs, DPPs and related objects, it seems standard to regard a refined 
enumeration as having a certain order, where this is based on only the number of boundary 
statistics being used. Hence, the enumeration of this paper is doubly-refined, although it 
involves two bulk statistics in addition to its two boundary statistics. 

The method of proving the main result of this paper is to introduce four-statistic, or 
doubly-refined, generating functions for ASMs and DPPs, and to show that these can be 
expressed as identical combinations of their respective three-statistic, or singly-refined, gen- 
erating functions. (The expression for the ASM case was previously obtained by Colomo 
and Pronko [ , Eq. (5.32)], [' ', Eq. (3.32)], but a new proof is given here.) The required 
equality of the ASM and DPP doubly-refined generating functions then follows from the 
equality, as confirmed in [0, Thm. 1], of the ASM and DPP singly-refined generating func- 
tions. The identities which give each doubly-refined generating function in terms of its 
singly-refined counterpart are obtained by applying the Desnanot-Jacobi determinant iden- 
tity to determinantal expressions for the generating functions. The latter expressions are 
themselves obtained using standard techniques involving the statistical mechanical six- vertex 
model with domain-wall boundary conditions and the Izergin-Korepin formula [31, 35] for 
the ASM case, and certain sets of nonintersecting lattice paths and the Lindstrom-Gessel- 
Viennot theorem [27, 28, l-'] for the DPP case. These techniques also played an essential 
role in the proof of [5, Thm. 1]. 
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An outline of the rest of this paper is as follows. In Section 2, the definitions of ASMs 
and DPPs, and of the relevant statistics and associated generating functions, are given. 
In Section 3, the main result involving the equality of the ASM and DPP doubly-refined 
generating functions, and the identities expressing each doubly-refined generating function 
in terms of its singly- refined counterpart, are stated. In Sections 5 and 6, proofs of the 
latter identity for the ASM and DPP cases respectively are given, with the Desnanot-Jacobi 
determinant identity, which is central to each of these proofs, having been introduced in 
Section 4. In Section 7, some further aspects of this work are discussed. 



2. Definitions 



2.1. ASMs and DPPs. In this section, the standard definitions of ASMs and DPPs, and 
of sets of ASMs and DPPs of order n, are given. 

An ASM, as first defined by Mills, Robbins and Rumsey [44, 45], is a square matrix in 
which each entry is 0, 1 or —1, the nonzero entries alternate in sign along each row and 
column, and the sum of entries in each row and column is 1. It follows that an ASM has a 
unique 1 in each of its first and last row and column, and that any permutation matrix is 
an ASM. 

A DPP, as first defined by Andrews [2, 3], is an array of positive integers, called parts, of 
the form 

DnDuDrs I^i,a, 

D22D23 D2 ,A2 + 1 

1^33 D, 



^3,A3+2 



(1) 



Dtt 



D 



t,Xt+t-l 



in which the parts decrease weakly along rows and decrease strictly down columns, and the 
first parts of the rows Da and row lengths Aj satisfy 



Dii > Ai > D22 > A2 > . . . > A-i,t-i > At_i > Du > Xf 

The empty array is also regarded as a DPP, and denoted 0. 

Examples of an ASM and a DPP (as also used in [H, Eqs. (1) & (4)]) are 



A 



/O 1 0\ 
1 0-110 
1-11 
1 
1 0-101 

VO 1 0/ 



6 6 6 5 2 
D= 44 1 . 
3 



(2) 



(3) 



For each positive integer n, the sets of ASMs and DPPs of order n are defined as 

ASM(n) = {nxn ASMs}, DPP(n) = {DPPs with each part < n}. (4) 
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For example, 

r/ioo\ /ooi\ /ioo\ /ooi\ /oio\ /oio\ /o i o\] 

ASM(3) = <^ 010, 010, 001,100,100, 001, l-lin>, 

[Vooiy VI 0/ Voioy Voioy Vooiy Viooy Vo i oyj 

DPP(3) = |0, ^ 2' 2, 3 3, 3, 3 2, 3 1 1 . (5) 



2.2. Statistics and generating functions for ASMs and DPPs. In this section, certain 
statistics, and associated generating functions, are defined for ASMs and DPPs. Some simple 
properties of the generating functions are also identified. 

For a given positive integer n, statistics for each A G ASM(?t,) are defined as 

^(^)= Yl ^v^^'f^ (6) 

l<i<i'<n 
i<j'<j<n 

fi{A) = number of — I's in A, (7) 
Pi{A) = number of O's to the left of the 1 in the first row of A, (8) 
P2{A) = number of O's to the right of the 1 in the last row of A, (9) 

and statistics for each D G DPP(r2) are defined as 

z/(D) = number of parts of D for which Dij > j — i, (10) 

p{D) = number of parts of D for which Dij < j — i, (11) 

Pi{D) = number of n's in D, (12) 

P2{D) = (number of (n — l)'s in D) + (number of rows of D of length n — 1). (13) 

It follows from (2) and (4) that, for the DPP statistic pi{D), n's can appear in only the first 
row of D, and that, for the DPP statistic p2{D), {n — l)'s can appear in only the first two 
rows of D, and that the number of rows of D of length n — 1 is or 1, since only the first 
row can have this length. 

For an ASM A, z/(y4) is regarded as the number of generalized inversions in A, since 
if A is a permutation matrix, i.e., if p{A) = 0, then z/(y4) is the number of inversions 
in the permutation vr given by 6-,^^^ = Aij. This statistic can also be written as z^(A) = 
Si<j<i'<n- i<j'<j<n ^i'j'^ whcrc this can be derived using the fact that the sum of entries 
in each row and column of A is constant. 

The definitions of the statistics u, p and pi for ASMs and DPPs are based on definitions of 
statistics introduced by Mills, Robbins and Rumsey [45, pp. 344-345]. There are some minor 
differences between the definitions used here and those of Mills, Robbins and Rumsey, with 
full details of these differences given in [5, Sec. 1.2]. For example, for an ASM or DPP X, 
Mills Robbins and Rumsey use i^{X) + p{X) instead of I'iX) as one of the basic statistics. It 
can be seen that, for an ASM A, i'{A)+p{A) can be regarded as an alternative generalization 
of the number of inversions in a permutation, since p{A) = if A is a permutation matrix. 
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For a DPP the parts Dij which contribute to ^{D) and /i(-D) in (lO)-(ll) are referred 
to as the nonspecial and special parts of D respectively, where this terminology matches that 
introduced by Mills, Robbins and Rumsey [ " , p. 344]. It can be seen that a DPP D contains 
at least as many nonspecial parts as rows, since the first part of each row is nonspecial. 

The statistic p2 for ASMs was first considered, in the context of its joint distribution 
with pi, and a conjectured connection with totally symmetric self-complementary plane 
partitions, by Mills, Robbins and Rumsey [ , p. 284]. The statistic p2 for DPPs has not 
previously been considered. 

It can be seen that, for an ASM or DPP X, v{X) and depend on the bulk structure 

of X, while Pi{X) and P2{X) depend only on the structure at or near the boundary of X. 

The statistics for the examples in (3) are, taking n = Q for the DPP D, 
u{A)=h, /i(A)=3, pi(A) = 3, P2{A) = 2 

u{D) = 7, piD) = 2, pi{D) = 3, P2{D) = 2. (14) 

Doubly-refined (or four-statistic) generating functions, which give weighted enumerations 
of the elements of ASM(r;,) or DPP(n) using arbitrary weights x, y, z\ and Z2 associated with 
the statistics (6)-(13), are now defined as 

zf'-(x,^,.„.2)= E (15) 

AeASM(n) 

Z°PP(x,i/,.„.2)= E (16) 

DGDPP(n) 

As indicated in Section 1, the term doubly-refined refers to the fact that these generating 
functions each involve two boundary statistics. 

As examples of (15)-(16), (5) gives 

Zf'^ix, 2/, zi, Z2) = Z3°PP(a;, z,, Z2) = 

1 + X^zlzl + XZ2 + x^z1z2 + XZi + x'^ziz'^ + XyZiZ2, (17) 
where the terms are written in an order which corresponds to that used in each set in (5). 
Singly- refined (or three-statistic) generating functions for ASMs and DPPs are now defined 

as 

Zf ^(x, y, z) = ^(x, y, z, 1), ^^^^(x, y, z) = ^^^^(x, z, 1). (18) 

The singly-refined and doubly-refined generating functions are also related, for n > 2, by 

Z^^'^ix.y.z) = Z„^s^(x,y,^,0), Z^ll{x,y,z) = ^^^^(x, y, z, 0). (19) 

These relations can be proved by constructing simple bijections between {A G ASM(r2 — 1) | 
v{A) = p, p{A) = m, pi{A) = k} and {A e ASM(n) | u{A) = p, p{A) = m, pi{A) = 
k, p2{A) = 0}, and between {D G DPP(n - 1) | u{D) = p, p{D) = m, pi{D) = k} and 
{D G DPP(n) I u{D) = p, p{D) = m, pi{D) = k, p2{D) = 0}. Specifically, an ASM in 
the first set is mapped to the second set by augmentation with a row and column on the 
bottom and right, where the common entry of the additional row and column is 1 and all 
other entries are 0, an ASM in the second set is mapped to the first set by deletion of its last 



6 



R. E. BEHREND, P. DI FRANCESCO, AND R ZINN- JUSTIN 



row and last column, a DPP in the first set is mapped to the second set by replacing each 
of its (n — l)'s by an n, and a DPP in the second set is mapped to the first set by replacing 
each of its n's by an n — 1. 

Similarly, it can be shown that Z^^Yi^^V^ ^) = Z^^^{x,y,0, z) and Z^^f{x,y,z) = 
y, 0, z), so that Z^^fix, 1) = Z^sm(^^ q, 0) and Z^^^x, y, 1) = Z^^^ix, y, 0, 0), 
although these relations are not needed for the main proofs of this paper. 



3. Main Results 

In this section, the main results of this paper are stated, and their connections to the main 
results of the paper [5] are pointed out. The overall method of proving the results is also 
outlined, but the presentation of the details of these proofs is deferred to later sections. 

It was conjectured by Mills, Robbins and Rumsey [ j , Conj. 3], and proved in [5], that 
the ASM and DPP singly-refined generating functions are equal. 

Theorem ([■'), Thm. 1]). For any n, x, y and z, 

Zf^(x,y,z) = ZDPP(x,i/,z). (20) 

A subsidiary result of Eqs. (97)-(98)] is that the singly-refined generating functions 
can be expressed as the determinant of a certain n x n matrix. Specifically, 

= = det {K^{x,y,z),,), (21) 

0<«j<n— 1 

where 

Kn[x,y, z)ii = —Oii+i + < , (22) 

" lEU„Et„(;:;.)("^i7>v-'^'. j=»-i- 

(For other, transformed versions of this determinant formula, based on formulae for DPPs 
obtained by Mills, Robbins and Rumsey [45, p. 346] and Lalonde [!0, Thm. 3.1], see, for 
example, [~, Eqs. (65)-(66) & (87)-(88)].) 

The primary result of this paper is that the ASM and DPP doubly-refined generating 
functions are also equal. 

Theorem 1. For any n, x, y, zi and z^, 

Zf^{x,y,z,,Z2) = Z°PP(a;,i/,zi,Z2). (23) 

Equivalently, for any n, p, m, ki and k2, the sizes of {A G ASM(n) | z/(A) = p, ^{A) = 
m, pi{A) = ki, p2{A) = k2} and {D e DPP(n) | u{D) = p, p{D) = m, pi{D) = 
ki, P2{D) = k2} are equal. 

It can be seen, using the definitions of ASMs and DPPs and the statistics (6)-(13), that 
each set in Theorem 1 is empty unless the integers p, m, ki and ^2 lie within the ranges 



(n-l)- 11 

p = 0,...,!^, m = 0,...,<'„(„l2)' k,,k2 = 0,...,n-l. (24) 



^ , n even. 
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A subsidiary result of this paper is that the doubly-refined generating functions can be 
expressed as the determinant of a certain n x n matrix. Specifically, 

Zn^^ix,y,zi,Z2) = Z^^^{x,y,zi,Z2) = _.det {Knix,y, zi, Z2)ij), (25) 

where 

Kn{x,y,Zi,Z2)ij = 

- + ELo Eto n-i>'y'-'^i^\ ^ = - - 2 (26) 

[eLoEIoELo Ct-,'Vf~'zT4-, J=n-1. 

(For an alternative, transformed version of this determinant formula see (65).) The matrices 
of (22) and (26) are related by Kn{x,y,z) = Kn{x,y,z, 1). 

Theorem 1 and (25) are valid for the trivial case n = 1 (for which ASM(l) = {(1)}, 
DPP(l) = {0} and Zf'^^^Xjy, Zi, Z2) = Z^^^{x,y, Zi, Z2) = 1). However, this case would 
require some distracting qualifications to be made to the statement of some subsequent 
results. It will therefore be assumed, throughout the remainder of this paper, that n > 2. 

Theorem 1 will be proved by deriving identities which express the ASM and DPP doubly- 
refined generating functions as identical combinations of their respective singly-refined gen- 
erating functions. The required equality (23) of the doubly- refined generating functions then 
follows immediately from the already-known equality (20) of the singly-refined generating 
functions. The identities expressing the doubly-refined generating functions in terms of their 
singly-refined counterparts are given by (27) in the following theorem. 

Theorem 2. Let Z denote Z^^^ or Z^^^ . Then, for any n, x, y, zi and Z2, 

{zi - Z2) Z„{x, y, zi, Z2) Zn-i{x, y, 1) = {zi - 1)2:2 Zn{x, y, zi) Z„_i(a;, y, Z2) - 

Zi{z2 - 1) Zn^i{x,y,zi) Zn{x,y,Z2). (27) 

Equivalently, for any n, x, y, Zi, Z2, z^ and z^, 

{Zi - Z2) (^3 - ^4) Zn{x, y, Zi, Z2) Znix, y, Z3, Z4) - 

{Zl - Z2,) {Z2 - Z4) Zn{x, y, Zi, Z3) Z„(x, y, Z2, Zi) + 

{Zi - Zi) {Z2 - Zs) Zn{x, y, Zi, Zi) Zn{x, ?/, Z2, Z3) = 0. (28) 

The ASM case of (27) with x = y = 1 was first obtained by Stroganov [-jS, Sec. 5], and 
the general ASM case of (27) was first obtained by Colomo and Pronko [12, Eq. (5.32)], [14, 
Eq. (3.32)]. In the latter case, (27) is expressed using the terminology of a two-point bound- 
ary correlation function for the six-vertex model with domain-wall boundary conditions. 
Related results for this boundary correlation function have also been obtained by Foda and 
Preston [25]. 

The equivalence of (27) and (28) can be verified as follows. Identity (28) follows from (27) 
by multiplying the LHS of (28) by Z„_i(x, y, 1)^, using (27) to express each of the six cases of 
{zi — Zj)Zn{x,y, Zi, Zj)Zn-i{x,y, 1) in terms of singly-refined generating functions, and then 
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checking that the resulting expression vanishes. Conversely, (27) follows from (28) by setting 
^3 = 1 and ^4 = 0, and then using (18)-(19). 

It also follows from (27) that 

Z^^^ix,y,z„Z2) = Z^'''{x,y,z,,z,), Z^^^ ix,y, z^, z,) = Z^^^ {x,y, z^, z,). (29) 

This symmetry in zi and Z2 can alternatively be derived directly using simple operations on 
ASMs and DPPs, as will be done in Section 7.2. 

Proofs of the ASM and DPP cases of Theorem 2 will be given in Sections 5 and 6 respec- 
tively. Each of these proofs will use the Desnanot-Jacobi determinant identity, which will be 
stated in Section 4. This identity will be applied to determinant formulae which will be ob- 
tained using standard methods, as also used in for parts of the proof of the equality (20) of 
the singly-refined generating functions. Specifically, for the ASM case of Theorem 2, a bijec- 
tion between ASMs and configurations of the six-vertex model with domain-wall boundary 
conditions, together with the Izergin-Korepin determinant formula for the partition func- 
tion of this model, will be used, and, for the DPP case of Theorem 2, a bijection between 
DPPs and certain sets of nonintersecting lattice paths, together with the Lindstrom-Gessel- 
Viennot theorem for the weighted enumeration of nonintersecting paths, will be used. The 
resulting determinant formula for the DPP case, which will be obtained in Section 6.3, will 
provide the second equality of (25). 

The proof given here of the ASM case of Theorem 2 will differ from the proof given by 
Colomo and Pronko [12, 11]. 

4. The Desnanot-Jacobi Identity 

In this section, the determinant identity, often referred to as the Desnanot-Jacobi identity, 
is stated, together with some background information. This identity will then be used in 
Sections 5 and 6. 

For a matrix M, and subsets / and J of the index sets for the rows and columns of M, 
let Mj denote the submatrix of M in which the rows of I and columns of J have been 
deleted. Omission of / or J in this notation means that the corresponding set is empty. 

Theorem (Desnanot-Jacobi). For any n x n matrix (Mjj) j^^. 1 < ii < i2 < n and 

1 < Ji < 32 < n, 

det M det M|;;;g = det M^J det M^} - det m|;^} det M^}. (30) 

Equivalently, for any {n+2)xn matrix {^ij) i^i^n+2-i< j<n ^^'^ < < ^2 < ^3 < ^4 < n+2, 
det A^i'^i''^^} jY{fe3,fe4} _ j\f{kiM} Ar{^2.'=4} ^ ^g^ ^{fci,fc4} ^^g^ j\j{k2M} = g. (31) 

In fact, only (31) will be directly used in this paper, but (30) is included here since it 
represents (usually with ii = ji = 1 and 12 = J2 = 2 or 12 = j2 = n) the more commonly- 
stated form of the theorem. 

The equivalence of (30) and (31) can be verified as follows. Identity (31) can be obtained 
from (30) by augmenting the (n + 2) x n matrix with two columns on the left, where 
these additional columns consist of I's in row ki of column 1 and row k2 of column 2 
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and O's elsewhere, and then applying (30) with ii = k^, 12 = k^, ji = 1 and j2 = 2 to the 
{n + 2) X (n + 2) matrix. Conversely, (30) can be obtained from (31) by augmenting the nxn 
matrix M with two rows on the top, where these additional rows consist of I's in column ji 
of row 1 and column j2 of row 2 and O's elsewhere, and then applying (31) with ki = 1, 
k2 = 2, ks = ii + 2 and k^ = 12 + 2 to the [n + 2) x n matrix. 

It can be shown similarly that (30) and (31) are also equivalent to the identity that, for 
any (n + 1) x n matrix (^ii)i<i<„+i.i<j<„, I < ki < k2 < k^ < n + 1 and I < I < n, 

det P^''^ det P/^/'"^^ - det P^^^> det Pg^''^^ + det P^'=3> det P/^^^'"^^ = 0. (32) 

An algebraic proof of (30) is given by Bressoud [8, Sec. 3.5], and a combinatorial proof of (30) 
is given by Zeilberger [62]. 

Cases of each of (30)-(32) for small values of n were published by Desnanot in 1819 (see 
Muir [47, Eqs. (A)-(G), (A')-(G'), pp. 139-142]). The further attribution to Jacobi is based 
on the fact that (30) for general n corresponds to the case m = 2 of the identity, published 
by Jacobi in 1834 (see Muir [17, Eq. (XX. 4), p. 208]), that for any nxn matrix M 
and any m < n, each size m minor of the matrix of size n — 1 minors of M equals the 
complementary minor of M multiplied by (det M)™"^. Proofs of the Jacobi identity, using a 
variety of methods, are given, for example, by Brualdi and Schneider [11, Sec. 4], Knuth [31, 
Eq. (3.16)], Leclerc [12, Sec. 3.2], Muir [18, Sec. 175] and TurnbuU [o'), pp. 77-79]. 

In fact, (30)-(32) also correspond to special cases of various other determinant identities. 
For example, (30) is a special case of an identity of Sylvester (see Brualdi and Berliner [10], 
Brualdi and Schneider [11, Eq. (8)], Knuth [31, Eq. (3.17)] or Leclerc [42, Eq. (8)]), and (31) 
is a special case of an identity of Bazin (see Leclerc [42, Eq. (7)]), which is itself a special 
case of an identity of Reiss and Picquet (see Leclerc [42, Sec 5.4]). Furthermore, (31) can be 
regarded as a simple case of a Pliicker relation (see, for example, Harris [29, pp. 65-66]). 

It is also interesting to note that the Desnanot-Jacobi identity is related to ASMs in a way 
seemingly-different to that established in this paper, namely through Dodgson's condensation 
algorithm for determinant evaluation [ ]. In this algorithm, the determinant of an n x n 
matrix M is computed by successively taking k = 2, . . . ,n, and using (30) (with ii = ji = 1 
and i2 = J2 = k) to compute all connected minors of M of size k, from the connected minors 
of M of sizes k — 1 and k — 2. If the coefficient of the second product of determinants on 
the RHS of (30) is changed to a parameter A, and this modified formula is instead used 
throughout Dodgson's algorithm, then the so-called A-determinant of M is produced, and 
remarkably, as shown by Robbins and Rumsey [ i )], this is naturally expressed as a sum over 
ASM(n), rather than as a sum over permutations of {1, . . . , n}. For further details see also, 
for example, Abeles [1], Bressoud [>^, Sec. 3.5] or Bressoud and Propp [9]. 

5. Proof of the ASM case of Theorem 2 

5.1. Bijection between ASMs and configurations of the six- vertex model with 

DWBC. In this section, the set of configurations of the statistical mechanical six- vertex 
model on an n X n grid with domain- wall boundary conditions (DWBC) is described, and the 
details of a natural bijection between ASM(n) and this set are summarized. This is standard 
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material, for which more information and references can be found in, for example, [j. Sees. 2.1 
& 3.1]. 

Let Qn be the nxn undirected grid with vertex set {(«, j) \ i, j = 0, . . . , ?t,+1}\{(0, 0), (0, n+ 
1), {n + 1, 0), [n + l,n + 1)}, where is taken to be in the ith row from the top and jth 
column from the left, and for which there are horizontal edges between {i,j) and ± 1), 
and vertical edges between (i, j) and {i ± 1, j), for each i,j = 1, . . .n. This grid is shown 
in Figure 1. The descriptions 'internal' and 'external' are applied in the obvious way to the 
vertices and edges of Qn- Specifically, the vertices of degree 4 are internal, the remaining 4n 
vertices of degree 1 are external, the edges which connect two internal vertices are internal, 
and the remaining edges are external. 



(1,0). 



(n,0)« 











(1,«) ^ 
































• (n,lj 
(n+ 


'.1) 




(n- 


(n,n) 
-l,n) 



«(l,ri + l) 



>(n,n + l) 



Figure 1. The undirected grid Qn- 



A configuration of the six-vertex model on Qn with DWBC is an assignment of arrows 
to the edges of Qn in which the external edge arrows on the upper, right, lower and left 
boundaries of Qn are all directed upward, leftward, downward and rightward, respectively, 
while the arrows on the four edges incident to any internal vertex satisfy the condition that 
two point towards and two point away from the vertex. 

Now define 

6VDW(n) = {configurations of the six-vertex model on Qn with DWBC}. (33) 
For example, 

6VDW(3) = 



(34) 



In an element of 6VDW(n), there are six possible configurations of arrows on the four 
edges incident to an internal vertex of Qn- These types of configurations, known as vertex 
configurations, are shown in Figure 2, where the numbers will be used to label the types. 

For any C G 6VDW(?t,), denote the total number of type-(/c) vertex configurations in C 
as Af(k){C), and denote the number of type-(A;) vertex configurations in C in row i of Qn 



as 
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(1) (2) (3) (4) (5) (6) 

Figure 2. The possible arrow configurations on edges incident to an internal vertex. 



It can be shown that, for any C G 6VDW(n), 

A/'(i)(C)=A/'(2)(C), A/-(3)(C)=Ar(4)(C), Ar(5)(C) 



A/'i)(C) 



A/'(5)(C) 



A/'(3)(C) 
1. 



A/'(6)(C) = 0, 



(35) 



Statistics for each C E 6VDW(n) are now defined as 

KC)=Ar(i)(C), /i(C^)=A/'(6)(C^), Pi(C)=Ar(\)(C), p,{C)=Ar^,^{C). (36) 

It can be shown straightforwardly that there is a natural bijection between ASM(ri) and 
6VDW(n), and that for each A G ASM(?2) and C G 6VDW(n) which correspond under 
this bijection, ^'{A) = z^(C), fi{A) = fi{C), pi{A) = pi{C) and P2(^) = P2{C) (where the 
statistics in these equations are given by (6)-(9) on the LHS and (36) on the RHS). 

The details of this bijection are as follows. To map A G ASM(?t,) to C G 6VDW(n), first 
associate the partial row sum Yl]'=i ^ij' with the horizontal edge between and + 1), 
for each i = 1, . . . ,n, j = 0, . . . ,n, and associate the partial column sum Yll'=i ^i'j with the 
vertical edge between and {i + 1, j), for each i = 0, . . . ,n, j = 1, . . . ,n. (The defining 
properties of ASMs imply that each of these partial sums is or 1.) Then obtain C by 
assigning a rightward or upward arrow to each edge associated with a 0, and a leftward or 
downward arrow to each edge associated with a 1. To map C G 6VDW(n) to A G ASM(n), 
set Aij to be 1, —1 or according to whether the vertex configuration in C at internal vertex 
(i, j) is of type (5), (6) or (l)-(4) respectively. 

For example, the ASM A of (3) and the corresponding six-vertex model configuration are 



/o 








1 


0\ 





1 





-1 


1 


1 


-1 


1 

















1 








1 





-1 


1 


Vo 








1 


0/ 



< — )■ 



(37) 



It can be seen that the statistics (36) for the six-vertex model configuration C in (37) are 
i^(C) = 5, /i(C) = 3, pi(C) = 3 and P2{C) = 2, and that these match the respective 
statistics (6)-(9) for the ASM A in (37), as given in (14). 

As further examples, in (5) and (34) the elements of ASM(3) and 6VDW(3) are listed in 
an order for which respective elements correspond under the bijection of this section. It can 
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be seen that Z^^^{x,y, Zi, Z2), as given in (17), could now be obtained using this bijection 
and (34). 

5.2. The partition function of the six-vertex model with DWBC. In this section, the 
partition function of the six-vertex model with DWBC is introduced. A relation between 
this partition function, at certain values of its parameters, and the doubly-refined ASM 
generating function (15) is then derived using the bijection of Section 5.1. 

Let a weight W(^k){u,v) be associated with each type (fc) of vertex configuration, where u 
and V are so-called spectral parameters. The partition function for the case of the six-vertex 
model of relevance here depends on these weights, and on parameters Ui and vj associated 
with row i and column j of Qn- Specifically, this partition function is defined as 

n 

C66VDW(n) i,j=l 

where (Cjj) is the type of vertex configuration in C at internal vertex 
Let the weights now satisfy 

W(i){u,v) = W(_2){u,v) = a{u,v), W(^3){u,v) = W(4.){u,v) = b{u,v), 

W(^5) {u,v) = W(^e){u,v) = c{u, v) , (39) 

for particular functions a{u,v), b{u,v) and c{u,v). 
If the spectral parameters in (38) are given by 

U2 = ... = M„_i = r, Ui= S, Un = t, Vi = . . . = Vn = W, (40) 

for parameters r, s, t and w, then 
Z^'^(s,r, . . . ,r,t;w, . . . ,w) = 



= 6(r, «;)("-i)("-2) ^) b(^t, w))"~^ c(r, w)"-^ c(s, w) c{t, w) x 

E/a(r,ui)\2A/'(i)(C) / c(r,to) \ 2A/'(6) (C) /a{s,w) b(r,to) \A^(\,(C) /a(t ,w) b{r,w) N •'Vjj) (C) 
\b{r,w) / \b{r,w)/ \ a{r,w) b{s,w) / \ a{r,w) b{t,w) / ' 

Ce6VDW(n) 

where (35) and the facts that Ylk=i-^{k)i^) — "^^ ^^"^ Ylk=i-^{k)i^) ~ were used. The 
bijection of Section 5.1, (6)-(9), (15) and (36) now give 

Z^V(s, r, . . . , r, t; w, . . . , w) = 6(r, w)("-i)("-2) (5(5^ ^) ^(t^ c(r, w)"^^ ^^^^ ^^^^ 

^ yASVlf / a(r,w) \2 / e(r,w) \ 2 a(g,'iD) b(7-,w) a(t,w) fc(r,w) \ / . ^ x 
n \Vf)(r,'U))/ ' U(r,w))/ ' a{r,-u;) 6(s,-!«) ' a(r,«)) fe{t,M;) T V / 
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5.3. The Izergin— Korepin determinant formula. In this section, the Izergin-Korepin 
formula for the partition function (38), with certain assignments of the weights (39), is 
stated. 

It was shown by Izergin [31], using certain results of Korepin [35], that if the weights (39) 
satisfy the Yang-Baxter equation, then the partition function (38) can be expressed in terms 
of the determinant of a certain n x n matrix. 

Let the weights (39) be given by 

a{u,v) = uq-^, b{u,v) = ^-vq, c{u,v) = {q"^ - j^) y/uv, (42) 

where g is a further, global parameter. These weights can be shown to satisfy the Yang- 
Baxter equation, with the resulting Izergin-Korepin determinant formula being given by the 
following result. 

Theorem (Izergin). The partition function (38) with weights (39) and (42) satisfies 

For a proof of (43) see, for example, Bogoliubov, Pronko and Zvonarev [7, Sec. 4], or 
Izergin, Coker and Korepin [32, Sec. 5]. 

It is apparent from (43) that Z^{ui, . . . , Vi, . . . , Vn) is symmetric in ui, . . . , Un and in 
f 1, . . . , f„. It can also be seen that, although the determinant and the denominator of the 
prefactor on the RHS of (43) both vanish if Ui = Uj or Vi = Vj for some i j, the RHS 
has a well-defined limit in these cases, since the LHS is well-defined for all such ui, . . . ,Un 
and vi, . . . ,Vn- Accordingly, it will be valid to use (43) to derive properties of the partition 
function for the values (40). 

5.4. Application of the Desnanot— Jacobi identity. In this section, the final steps in the 
proof of the ASM case of (28) are taken. These involve using the form (31) of the Desnanot- 
Jacobi identity, the Izergin-Korepin determinant formula (43), and the relation (41) between 
the doubly-refined ASM generating function and the partition function of the six-vertex 
model with DWBC. 

Applying the form (31) of Desnanot- Jacobi identity to the [n + 2) x n matrix given 
by \ -t—^'-tt/^ — \j ; and then applying the Izergin-Korepin determinant for- 

\a(u^,Vj) t)(Ui,Vj)j ]^<j<„_,_2. l<j<n 

mula (43) to each of the six determinants which appear, it follows that, for any ui, . . . , Un+2, 
Vi, . . . ,Vn and 1 < ki < k2 < < < n + 2, 

{Uki - Mfca) («fc3 - ^^4) ^n^(^l5 • • • . . .,Uk2, ■ ■ ■,Un+2;Vl, ■■■,Vn) X 

Z^^{ui, . . . , Mfc3, . . . , Ufc4, . . . , Un+2; f 1, • • • , Vn) 
- {Uki - Mfca) (Wfca - Uk^) {ui, . . . , Mfc^ , . . . , , . . . , Un+2] f 1, • • • , ^n) X 

Zf^iUi, . . . , Mfc2, . . . , . . . , Un+2] f 1, • • • , Vn) 

+ {Uki - UkJ {Uk2 - Uks) Z^^{ui, . . .,Uki, ■ ■ ■,Uki, ■ ■ ■,Un+2; Vl, ...,Vn) X 

Zn^iUi, . . . , . . . , Mfe3, . . . , Un+2; Vl,..., Vn) = 0, (44) 
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where Ui, . . . ,Ui, . . . ,Uj, . . . , Un+2 denotes the omission of Ui and Uj from ui, . . . , Un+2- Note 
that (44) would be satisfied by any function which has the form of (43), for arbitrary functions 

a{u,v), b{u,v) and c{u,v). 

Setting M5 = . . . = Un+2 = r, vi = . . . = Vn = w and fcj = i in (44) gives 

{ui - U2) (m3 - M4) Z^{ui, r, . . . ,r,U2;w, . . . ,w) {u^., r, . . . , r, M4; w, . . . , w) - 

{ui - M3) {U2 - Ui) Z^{ui, r, . . . ,r,U3;w, . . . ,w) 2'^^(m2, r, . . . , r, M4; w, . . . , w) + 
{ui - U4) {U2 - U3) Zf^{ui, r,...,r,U4;w,...,w) Z^^ {u2, r, . . . , r, U3; w, . . . , u;) = 0, (45) 

where the symmetry of {u\, . . . , Un\ vi, . . . , fn) in mi, . . . , m„ has been used. 
Now let parameters x, y and zi, . . . , Z4 he given, using the functions of (39), as 

_ / a(r,«)) \2 _ / c(r,w) \2 _ a(n^,-»)) fc(r,«)) • _ . /.^x 

It can be checked that this allows arbitrary x, y and zi, . . . , Z4 to be expressed in terms of 
parameters q, r, w and Ui, . . . , U4. (In fact, this remains the case for rw = 1.) It can also be 
checked that, for z, j = 1, . . . , 4, 

„ - a-{r,w)b{ui,w)b{uj,w) / x /.„x 

~ 6(r,W))(g-2-g2)«, ^ ' 

The ASM case of (28) now follows from (45), using (41), (47) and the fact that the arbi- 
trary X, y and zi, . . . , ^4 in (28) can be parameterized by (46). 

It is interesting to note that the Desnanot-Jacobi identity was also used together with 
an Izergin-Korepin determinant formula by Korepin and Zinn- Justin [ , Sec. 3] (see also 
Sogo ["'>, Sec. 4]) in the computation of the thermodynamic limit of the six- vertex model 
with DWBC in its ferroelectric and disordered regimes. 



6. Proof of the DPP case of Theorem 2 

6.1. Bijection between DPPs and sets of nonintersecting lattice paths. In this 
section, the details of a bijection between DPP(n) and a set of certain sets of nonintersecting 
lattice paths on an n x grid are summarized. More information, background material and 
references can be found in [ , Sees. 2.2 & 3.2]. 

Let Qn be the nxn directed grid with vertex set \ i, j = 0, . . . , n—1}, where is in 

the ith column from the left and jth row from the bottom, and for which there are horizontal 
edges directed rightward from to (i + 1, j) for i = 0, . . . , — 2, j = 0, . . . , — 1, and 
vertical edges directed downward from (i, j) to (z, j — 1) for i = 0, . . . , n — 1, j = 1, . . . , n — 1. 
This grid is shown in Figure 3, where the edge weights in this diagram will be introduced 
in Section 6.3. Note that, in contrast to Qn in Section 5.1, the rows and columns of Qn are 
labeled by Cartesian, rather than matrix-type, coordinates. 
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Figure 3. The directed grid with associated edge weights. 



Now define 
NILPfn) = 



sets P of noninter- 
secting paths on Qn 



there exist < t < — 1 and n = \q > \i > . . . > 
Xt > Xt+i = for which P consists of paths from 
(0, A,;_i — 1) to (Aj, 0), for each i = 1, . . . ,t + 1 



For example, 

NILP(3) 



Statistics for each P G NILP(r;,) are now defined as 

z/(P) = number of rightward steps of P above the subdiagonal fine — 1)}, 
yu(P) = number of rightward steps of P below the subdiagonal line — 1)}, 
Pi{P) = number of rightward steps of P in row n — 1 oi Qn, 
P2{P) = (number of rightward steps of P in row n — 2 of Qn) + 
(number of paths of P which start at (0, n — 2)). 



(4J 



(49) 



(50) 
(51) 
(52) 

(53) 



Note that the second term on the RHS of (53) is also the number of paths of P which end 
at (n — 1, 0). 

It can be shown straightforwardly that there is a natural bijection between DPP(n) and 
NILP(n), and that for each D G ASM(?t,) and P G NILP(r;,) which correspond under this 
bijection, u{D) = u{P), p{D) = fi{P), pi{D) = pi(P) and p2{D) = p2{P) (where the 
statistics in these equations are given by (10)-(13) on the LHS and (50)-(53) on the RHS). 

The details of this bijection are as follows. To map D G DPP(n) to P G NILP(n), first 
let t be the number of rows in D and let Aj be the length of row i of as in (1). Also define 
\q = n and \t+i = 0. Then obtain P by forming a path, for each i = 1, . . . ,t + 1, from 
(0, Ai_i — 1) to (Aj, 0) whose kth rightward step has height Dj^j+fc-i — 1- To map P G NILP(n) 
to D G DPP(n), set Dij to be 1 plus the height of the {j — i + l)th rightward step in the ith 
path from the top. 
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For example, the DPP D of (3) (with n = 6) and the corresponding set of nonintersecting 
lattice paths are 



66652 
44 1 
3 



i — )■ 




(54) 



where the subdiagonal line — 1)} is shown, and each part of D is displayed above 
its corresponding rightward step, with the nonspecial and special parts in red and green 
respectively. It can be seen that the statistics (50)-(53) for the set of paths P in (54) are 
i/(P) = 7, /i(-P) = 2, pi(-P) = 3 and P2{P) = 2, and that these match the respective 
statistics (10)-(13) for the DPP D in (54), as given in (14). 

As further examples, in (5) and (49) the elements of DPP (3) and NILP(3) are listed in 
an order for which respective elements correspond under the bijection of this section. It can 
be seen that Z^^^{x,y, zi, Z2), as given in (17), could now be obtained using this bijection 
and (49). 



6.2. The Lindstrom— Gessel— Viennot theorem. In this section, the Lindstrom-Gessel- 
Viennot theorem for the weighted enumeration of sets of nonintersecting paths in terms of 
a determinant is stated. 

Consider an acyclic directed graph G, let a weight be assigned to each edge of G, and 
define the weight W{p) of a path p on G to be the product of the weights of the edges along 
which p passes. For vertices u and v of G, let Vu,v denote the set of all paths on G from u to v. 
For vertices Ui, . . . , Um and f 1, . . . , of G, let Afdui, . . . , Um', f 1, . . . , Vm) denote the set of 
all sets P of paths on G such that P consists of a path of Vu,,vi for each i = 1, . . . , m, and 
different paths of P are nonintersecting. The Lindstrom-Gessel- Viennot theorem [27, 28, 43] 
can now be stated as follows. 

Theorem (Lindstrom; Gessel, Viennot). // Afciuaj^, . . . , Ua^; vi, . . . , Vm) is empty for each 
permutation a of {1, ... , m} other than the identity, then 

For a proof of (55) see, for example, Gessel and Viennot [-u. Sec. 2] or Stembridge [57, 
Sec. 1]. 

6.3. Determinant formula for the DPP doubly-refined generating function. In this 
section, the bijection of Section 6.1 and the Lindstrom-Gessel- Viennot theorem (55) are used 
to obtain a determinant formula, as already stated in the second equality of (25), for the 
doubly-refined DPP generating function (16). 
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The directed graph G of Section 6.2 is now taken to be the n x directed grid defined 
in Section 6.1. 

The set NILP(n), as defined in (48), can be written using the notation of Section 6.2 as 
NILP(n) = 

U Mg^ ((0, n - 1), (0, Ai - 1), . . . , (0, Ai - 1); (Ai, 0) . . . , (A*, 0), (0, 0)). (56) 

0<t<n-l 
n-l>Ai>...>At>l 

Now assign weights to the edges of as indicated in Figure 3. More specifically, assign xzi 
to each horizontal edge in row n — 1, xz2 to each horizontal edge in row n — 2, x to each 
horizontal edge from («, j) to {i + 1, j) for < i < j < n — 3, ?/ to each horizontal edge from 
(z, j) to {i + 1, j) for < j < i < n — 2, and 1 to each vertical edge. Also, shghtly expanding 
the notation of Section 6.2, denote the weight of a path p on Qn as W{x, y, zi, Z2,p). 

The bijection of Section 6.1 and the Lindstrom-Gessel-Viennot theorem (55), together 
with (16), (50)-(53), (56) and the observation that the condition for the validity of (55) is 
satisfied, now give 

Zn^^{x,y,zi,Z2) 

{x,y,zi,Z2,p) 1 I . (57) 




det ( 

i=At-l,...,Ai-l,n-l ^ Pe^(Oj).(>,0) 



Note that the term Z2^^ ^ (= -^2^^^ " ^) ^^e RHS of (57) arises from the second term 
on the RHS of (53). 

It can be shown straightforwardly that, for any matrix (Mjj)^^. 

^^det^_^(Mij - = ^ detM{o}us,(5--i)u{n-i}, (58) 

S'C{l,...,n-l} 

where M{o}us',{s-i)u{n-i} denotes the submatrix of M formed by restricting the rows and 
columns to those indexed by {0} U S and {s — 1 | s G S*} U {n — 1} respectively. 

Applying (58) to (57), and taking account of the term ^2'^^^"^'"~^ gives 



^ W{x,y,zi,Z2,p)\ 

^(0,j),(i,0) / 



Z^^^{x,y,Zi,Z2)= ^det + >J W{x,y, z,, Z2,p) ] . (59) 

(Note that Z^^^{x, y, zi, Z2) = deto<ij<„-i(-(5»j+i + Z2*-"~' T.per^^,^,^,^,^o)^i^^ ^i' ^2,p)) is 
also valid.) 
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j < n — 3 
j = n-2 
'4"'", J = n-1. 



(60) 



It is now found that 
^ W{x,y,Zi,Z2,p) = 

Z]fc=o Z]i=o (i-fc) ( k-i )^^y'' ''4> 

sr^k sr^l n-i\,n-l-2\ ki-km 
. Z^fc=0 Z^«=0 Z^m=0 U-fcA k-l /-^ y ^1 

More specifically, the sums of weights of paths of 'P(o,j),{i,o) in (60) can be obtained as follows, 
with these derivations being shown diagrammatically in Figure 4. For j < n — 3, combine 
any of the {^^) paths of P(oj),(A;,fc-i), each with weight x^, and any of the (*~^) paths of 
(j o), each with weight y^~^., for any < < min(i, j + 1). For j = n — 2, combine 
the single path of 'P(o,n-2),(/,n-2), with weight (^2:2)', the single path of 'P(i,n-2),(z,n-3), with 
weight 1, any of the ('^^17^) P^-^hs of 'P(z,n-3),(A:,fc-i), each with weight x''~\ and any of the (*~^) 

1, combine 



n 



paths of V(k,k-i),{i,o), each with weight ?/* , for any < / < A; < z. For j 
the single path of P(o,n-i),(m,n-i), with weight [xzi)"^, the single path of V{m,n-i),{m,n-2), 
with weight 1, the single path of 'P(m,n-2),(;,n-2), with weight (xz2)'~™, the single path of 
^(«,n-2),o,n-3), with Weight 1, any of the ("^i7^) paths of P(/,n-3),(fc,fc-i), each with weight x''~\ 
and any of the (Izl) paths of 'P(fe,fc-i),(j,o); each with weight y^~'^, for any 0<m<l<k<i. 




n-l: 



n-1 




n-l I 




n—l m I k i n-l 



j < n — 3 



J =n-2 
Figure 4. Derivation of (60). 

Finally, using (60) in (59) gives the second equality of (25), i.e., 

Zn^^ix,y,zi,Z2) =_.det (Kn{x,y,zi,Z2)ij), 

0<«j<n— 1 

where Kn{x,y, zi, Z2)ij is defined in (26). 



J = n 



(61) 



6.4. Transformation of the DPP determinant formula. In this section, a further de- 
terminant formula for the DPP doubly-refined generating function (16) is obtained by ele- 
mentary transformation of the determinant formula (61). 

Define, for i = 0, . . . ,n — 1, 



c„(x, y, Z2). = ELo Eto ELo (::D rk-'>'y'-'^i 



J—m 



(62) 
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Thus, Cn{x, y, Z2)i — 5i,n-i and y, Zi, ^2)1 are the entries in row i of the second-last and 

last columns respectively of the matrix Kn{x, y, Zi, Z2), as defined in (26) and used in (61). 

It is now found that 

{zi - Z2)Cn{x, y, zi, Z2)i = C„(x, y, zi)i - Cn{x, y, Z2)i. (63) 
This can be proved combinatorially as follows. As shown in Section 6.3, and stated in the 
cases j = n - 1 and j = n - 2 of (60), C„(x, zi, Z2)i = EpGP(o,„_i),(i,o)^(^' ^' ^i' ^2,p) and 
Z2~'^ Cnix,y,Z2)i = X]pe7'(o,„-2).(>,o)^(^'^'^i'^2,p)- By partitioning P(o,n-i),(i,o) into those 
paths with an initial downward step, which correspond (by deleting the initial step, with 
weight 1) to 'P(o,n-2),(i,o)) and those paths with an initial rightward step, which correspond 
(by deleting the initial step, with weight xzi) to 'P(i,n-i),(j,o), it follows that 

{xzi)~'^{Cnix,y,zi,Z2)i- Z2~'^ Cn{x,y,Z2)i) = ^ W{x,y,zi,Z2,p). (64) 

Now consider the involutions on 'P(o,n-i),(i,o) and on P(i,n-i),(j,o) in which a path which starts 
with ki rightward steps in row n — 1 of followed by ^2 rightward steps in row n — 2 
of Qn, is mapped to the path in which the roles of ki and ^2 are interchanged, with the part 
of the path below row n — 2 being kept the same. Using these involutions, it follows that 
Cn{x,y, zi, Z2) and the LHS of (64) are symmetric in zi and Z2, which gives (63). 

Alternatively, (63) can be proved algebraically as follows. Let [m*]/(m) denote the coeffi- 
cient of in a power series f{u). Then, briefly, C„(x, y, z)i = [u'] jz^^^^^i^ + T^)"~ > 
and Cn{x, y, z,), - Cn{x, y, Z2)^ = [u^] ( - J (1 - yu) (l + i^)""' = 

{l-ixizi-l)+y)l)il-ixiz2-l)+y)u) + ^ (^1 ~ Z2)Cn{x, y, Zi, Z2)i- 

Starting with the determinant formula (61), together with definitions (26) and (62), and 
then successively multiplying the last column of K^^x, y, Zi, Z2) by 21 — 2:2, adding the second- 
last column to the last column, applying (63), and interchanging the last two columns, it 
now follows that 



{z2-zi)Z^^^{x,y,zi,Z2)= dei (L„(x, y, 21, 22)*^), (65) 

0<«,.j<ra— 1 

where 



0<ij<r 

eS^'^'^'^ (::;) tuVy'-' - J < - - 3 

:„(x,y,2i,22)ij = C„(x,?/,2i)i - j = n-2 (66) 

Cn{x,y,Z2)i - 6i^n-l, j = H - 1. 

6.5. Application of the Desnanot— Jacobi identity. In this section, the final steps in 
the proof of the DPP case of (28) are taken. 

Define an n x (n + 2) matrix (L„(x, y, 21, 22, 23, Z4)ij) Q^i^^_^.Q^j^^^^ by 

{Ln{x,y,zi,Z2)ij, j<n-l 
Cn{x, y, 23)i - 6i^n-i, 3 = n (67) 
C„(x, 2/,24)i - 5i,„_i, j = n + l, 

where Ln{x,y, zi, Z2)ij and Cn{x,y, z)i are given by (66) and (62). 
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Applying the form (31) of the Desnanot-Jacobi identity with its four selected rows taken 
aski = n + i — 3 (and with row and column numbers now starting from 0) to the transpose 
of the matrix (67), and using the determinant formula (65) on each of the six determinants 
which appear, now gives the DPP case of (28). 

7. Discussion 

The proof of Theorem 1, which has been the primary focus of this paper, has now been 
presented. In summary, the proof involved using certain bijections together with the Izergin- 
Korepin formula (43) and Lindstrom-Gessel-Viennot theorem (55) to obtain determinan- 
tal expressions for the doubly-refined enumeration of ASMs and DPPs, then applying the 
form (31) of the Desnanot-Jacobi identity to obtain the identity (28) satisfied by the ASM 
and DPP doubly-refined generating functions, then using the argument given in Section 3 
(i.e., setting 2:3 = 1 and = in (28) and applying (18)-(19)) to show that the doubly- 
refined generating functions can be expressed in terms of their singly-refined counterparts 
by (27), and finally using the equality (20) of the singly- refined generating functions to 
establish the required equality (23) of the doubly-refined generating functions. 

A variation of this proof involves applying the form (32), instead of (31), of the Desnanot- 
Jacobi identity to the ASM and DPP determinants, which leads directly to (27), instead 
of (28). However, in this approach, different terms of (27) need to be treated separately 
at certain stages of the derivation, whereas the approach which was used seems preferable, 
since all of the terms of (28) can be treated uniformly throughout. 

In the remainder of this section, some further matters related to the main results of this 
paper are discussed. 

7.1. Special cases. In this section, some special cases of Theorem 1 which can be proved 
bijectively are discussed briefly. 

It follows from Theorem 1 that, for each n, there exist bijections between ASM(n) and 
DPP(n) with the property that, for each A G ASM(n) and D G DPP(r;,) which correspond 
under the bijection, z/(y4) = z^(-D), Ai(v4) = fi{D), pi{A) = pi(-D) and p2{A) = p2{D). 

However, although these bijections necessarily exist, a natural and explicit such bijection 
for arbitrary n is not currently known. Nevertheless, explicit bijections are known for some 
special classes of the subsets of ASM(n) and DPP(n) in Theorem 1, i.e., for certain classes of 
values of p, m, ki and /c2 in Theorem 1, with these bijections thereby providing alternative 
proofs of certain cases of the theorem. Some of these cases will now be outlined, although 
most of the details will be omitted. 

For ki = k2, p = ki{ki + l)/2 and m = ki{n — fci — 1), the sets in Theorem 1 consist of a 
single ASM and single DPP, for each 0<A;i = A;2<'^ — 1 and arbitrary n, so that there is 
a single, trivial bijection. The ASM and DPP are given explicitly in [ j, Eqs. (21)-(22)]. 

For m = 0, with n, p, ki and ^2 arbitrary (i.e., for permutation matrices and DPPs 
with no special parts), a bijection can be obtained by modifying a known bijection for the 
corresponding singly-refined case in which p2 and /c2 do not appear, where the modification 
is needed to give the property that P2{A.) = P2{D) (= k2) for any corresponding ASM A and 
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DPP D. Details and references for the known bijection are given in [ ), pp. 336-337]. It can 
also be shown straightforwardly for this case that 

^(x, 0, ^1, z,) = 0, ^1, z,) = 

[n-2]J {x^'+'-^-^ziz^-^-^ + x''-'+'-h'^-'^^zi), (68) 

0<i<i<n-l 

where [n]^. = 1 + x + . . . + x""^ and [n]J = [n]r,[n - 1]^. . . . [l]^;. 

Forp = 1, with n, m, ki and ^2 arbitrary, a bijection can again be obtained by modifying 
a known bijection for the corresponding singly-refined case in which p2 and k2 do not appear, 
with the modification needed to ensure that p2{A) = P2{D) for any corresponding ASM A 
and DPP D. Details for the known bijection are essentially given in [5, p. 337]. This case, 
together with the operation * which will be described in Section 7.2, can also be used to give 
a bijection for the case p = n{n — l)/2 — m — 1, with n, m, ki and ^2 arbitrary. 

In these cases, the bijections for the singly-refined cases seem reasonably natural, whereas 
the subsequent modifications do not. In fact, it seems that a natural, full bijection between 
ASM(n) and DPP(n) might satisfy u{A) = u{D), fi{A) = p{D), pi{A) = pi{D) and p2{A) = 
P2{(t{D)) for each corresponding A G ASM(n) and D G DPP(n), where cr is a permutation 
(other than the identity) of DPP(n) which satisfies z/((t(D)) = i^iD), p{a{D)) = p{D) and 
Pi{a{D)) = pi{D) for each D G DPP(n). If this is so, then it would be natural to replace 
the fourth statistic for each D G DPP(r2) by p2{(y{D)). 

7.2. Symmetry operations. In this section, two symmetry operations on ASM(n) and 
DPP(n) which lead directly to some common properties of the doubly- refined generating 
functions Z^^^{x,y, Zi, Z2) and Z^^^ {x,y, zi, Z2) are discussed. 

Operations * and f which transform A G ASM(?2) to A*, A''' G ASM(n) are defined by 

= ^i.n+l-j; = ^n+l-i,n+l-j- (69) 

Thus, * and f correspond to reflection of A in a vertical line, and rotation of A by tt 
respectively. These operations, together with the other natural reflections and rotations 
of ASMs, were first considered by Mills, Robbins and Rumsey [15], [46], [53], [54], who 
formulated numerous conjectures involving the sizes of sets of ASMs invariant under groups 
of such operations. Many of these conjectures have since been proved, and further related 
results have been obtained, by Kuperberg [■)*)], Okada [19], and Razumov and Stroganov [51, 
52]. 

An operation * which transforms D G DPP(n) to D* G DPP(n) is defined as follows. 
If Dn-j,n-i is defined and Dn-j,n-i < j — then D*^ = j — i + 1 — Dn-j,n-i, if i < j < n 
and Dn~j,n-i is not defined then D*j = n + 1 — i — (number of positive integers k satisfying 
n + 2 — z — Dk^n-j < k < n — j), and if Dn-j^n-i > 3 — i then D*- is not defined. This 
operation was first defined by Mills, Robbins and Rumsey [45, p. 351], using the previous 
rather complicated description. However, it was shown by Lalonde [ ! 1] that the operation 
has a much simpler description in terms of the set of nonintersecting lattice paths which 
correspond to a DPP. (Lalonde used slightly different nonintersecting paths from the paths 
of (48) used here, but it can be shown straightforwardly that * also has a relatively simple 
description in terms of (48).) Furthermore, it was shown by Krattenthaler [37] that there 
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is a natural bijection between DPP(n) and the set of cyclically symmetric tilings with three 
types of unit rhombi of a certain punctured hexagon with alternating sides of length n — 1 
and n + 1, and that, in terms of such tilings, * simply corresponds to the reflection of a tiling 
in a symmetry axis of the hexagon. 

An operation f which transforms D G DPP(?t,) to D'^ G DPP(n) is defined as follows. 
Let Pi(-D) and P2{D) be given by (12)-(13) and pz{D) be the number of {n — l)'s in the first 
row of D. Then obtain D'^ by replacing the Pi(-D) ^^-'s and Pz{D) {n — l)'s in the first row 
of D by P2{D) ri's followed by (pi(-D) + PsiD) — P2{D)) {n — l)'s, leaving all other parts of D 
unchanged. This operation has not previously been considered in the literature. 

It can be seen that the operations * and f are involutions on ASM(n) and on DPP(ri). 
Furthermore, the statistics (7)-(13) behave under these operations according to 

^(X*) = ^ - u{X) - p{X), p{X*) = p{X), 

pi{X*) = n - 1 - pi(X), p2(X*) = n - 1 - p2(X), (70) 

and 

u{X^) = u{X), p(Xt)=p(X), pi(Xt)=p2(X), (71) 

for any X G ASM(?2) or X G DPP(n). Proofs of the first three cases of (70) for ASMs 
and DPPs are given by Mills, Robbins and Rumsey [ 10, p. 352], and for DPPs by Kratten- 
thaler [37, pp. 1143-1144] and Lalonde [11, p. 317]. The remaining cases in (70)-(71) can 
be obtained similarly. 

It follows from (70) that *, when restricted to the set of ASMs or the set of DPPs in 
Theorem 1, provides a bijection to the corresponding set in which p, ki and /c2 are replaced 
by n(n — 1) /2 — p — m, n — 1 — ki and n — 1 — A;2 respectively. Therefore, 

^(x, y, z„ Z2) = {z.z^r-' ^(i, f , ^, ^), 

Z°PP(x, y, z,, Z2) = {z.z^r-' Z°PP(i, y-, ^, ^). (72) 

Similarly, it follows from (71) that f, when restricted to the set of ASMs or the set of 
DPPs in Theorem 1, provides a bijection to the corresponding set in which ki and k2 are 
interchanged. Therefore, Z:^^^{x, y, Zi, Z2) and Z^^^{x, y, zi, Z2) are symmetric in zi and Z2, 
as already noted in (29). 

The question of whether * and f for ASMs can correspond to * and f for DPPs under a 
bijection between ASM(n) and DPP(n) will now be addressed. Since * and f are involutions, 
whose orbits are therefore of size 1 or 2, this question reduces to a question of whether the 
numbers of ^-invariant or f-invariant objects are equal for ASM(ra) and for DPP(n). 

In the case of *, it follows from a result of de Gier, Pyatov and Zinn- Justin [ , Prop. 3, first 
equation], initially conjectured by Mills, Robbins and Rumsey [15, Conj. 3S], that the sizes of 
the *-invariant sets {A G ASM(n) \ A* = A, v{A) = p, p{A) = m, pi{A) = ki, p2{A) = ^2} 
and {D G DPP(n) \ D* = D, v{D) = p, p{D) = m, pi{D) = ki, p2{D) = ^2} are equal 
for any ri, m, ki and /c2- (It can be seen, using (70) and the structure of *-invariant 
ASMs and DPPs, that each of these sets is empty unless n is odd, 2p + m = n{n — l)/2, 
^1 = ^2 = (^ — l)/2, and m/2 — (n — l)/4 is a nonnegative integer.) It now follows from 
this result and Theorem 1 that there exist bijections between ASM(n) and DPP(n) with 
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the properties that, for each A G ASM(n) and D G DPP(n) which correspond under the 
bijection, A* and D* also correspond under the bijection, z/(y4) = z^(-D), fJ'iA) = fJ,{D), 
pM) = Pi{D) and p^) = p^^D). 

In the case of f, the last equation of (71) implies that if an ASM or DPP X is f-invariant 
then pi{X) = p2{X). Furthermore, it can be seen using the definitions of f that, for n > 4, 
there exist non- f-invariant n x n ASMs A which satisfy pi{A) = P2{A), but that for all n, 
D G DPP(n) is f-invariant if and only if Pi(-D) = P2{D). Therefore, for n > 4 there does not 
exist a bijection between ASM(n) and DPP(?2) with the property that, for each A G ASM(n) 
and D G DPP(n) which correspond under the bijection, A^ and also correspond under 
the bijection. However, this is not surprising, since f for an ASM involves its bulk structure, 
whereas f for a DPP involves only its top two rows. 

7.3. Further multiply-refined enumeration results. In this section, other enumerative 
results for ASMs, and in some cases DPPs, involving one or more boundary statistics, but 
mostly no bulk statistics, are discussed briefly. 

An ASM contains four natural boundary statistics corresponding to the positions of the I's 
in its first and last row and column. More specifically, for an ASM A, these statistics can be 
taken as pi{A), . . . , p4{A), where pi{A) and P2{A) are given by (8)-(9), and 

PsiA) = number of O's above the 1 in the first column of A, 

p4{A) = number of O's below the 1 in the last column of A. (73) 

It can easily be seen, using the natural symmetry operations on ASMs, that the joint 
distribution on ASM(?t,) of the bulk statistics u and p (as given by (6)-(7)) together with 
a single boundary statistic pi is independent of i, for i = 1,...,4. It follows from the 
result (20), together with (18) and (29), that these distributions on ASM(n) also equal 
the joint distribution on DPP(n) of u, p and pj (as given by (10)-(13)), for j = 1,2. A 
determinant formula for these six identical distributions is provided by (21). 

Now consider singly-refined enumeration without bulk statistics, and define 

An,k = \{A G ASM(n) I p,{A) = k}\ = \{D G DPP(n) | p,(D) = k}l (74) 
for any 1 < i < 4 and 1 < j < 2. Then 

A _ (n+fc-l)! (2n-fc-2)! ■r-rn-2 (31+1)] 
-Hn.fc — (2„„2)! fc! (n-fc-l)! 1 h=0 (n+i-1)! ' ^'"^^ 

where this formula for the DPP case follows from results of Mills, Robbins and Rumsey [ t l] 
(see Bressoud [ , Conj. 9 and Sec. 5.3] for further details), while the formula for the ASM case 
was first proved by Zeilberger [Gl] (see Bressoud [ , Sec. 7.3] for further details), following 
conjectures of Mills, Robbins and Rumsey [11, Conj. 1], [15, Conj. 2]. Alternative proofs of 
the ASM case have been given by Colomo and Pronko [ , Sec. 5.3], [11, Sec. 4.2], Fischer [_u] 
and Stroganov [oN, Sec. 4]. Each of these proofs of the ASM case of (75) uses the six- vertex 
model with DWBC and particular forms of the Izergin-Korepin formula, except for the 
proof of Fischer [20], which makes essential use of a certain operator formula obtained by 
Fischer [19, 21]. Now define An = |ASM(n)| = |DPP(n)| = Efc^ Then, as discussed 
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in Section 1, and as can be obtained from (75), 

A = TT""^ (^'+^' (76) 

Proceeding to doubly-refined enumeration, it can easily be seen that the joint distributions 
on ASM(ra) of the four statistics z/, fi, pi and p2, and the four statistics fi, ps and p^ (as given 
by (6)-(9) and (73)) are equal (and that each is symmetric in its two boundary statistics), 
and Theorem 1 states that these distributions are also equal to the joint distribution on 
DPP(n) of u, p, pi and p2 (as given by (10)-(13)). A determinant formula for these three 
identical distributions is provided by (25). 

Now consider doubly-refined enumeration without bulk statistics, and define 

An,^,J = \{Ae ASM(n) I pi(A) = t, p2{A) =j}\ = \{Ae ASM(n) | p^{A) = t, p^{A) = j}\ 

= \{De DPP(n) I pi(D) = 2, p2{D) =j}\. (77) 

It follows from (27) that 

('^riji— IJ •^n,i,j — l) ^n~l ^n,j— 1 ^n— — 1 ^n,i '^n— 1 j'— 1 

•^n— 1 •A-nJ — l "^rt— 1 •A-nj, C^S) 

which can be solved to give 
1 



^ min(i,n— jr — 1) 



'^n— fc — 1 '^nj'+fc+l ~l~ '^n— fc— 1 j'+fe) • ('''9) 

Thus, (79), together with (75)-(76), provides an explicit formula for ^n,ij- The ASM case 
of (78) was first obtained by Stroganov ['s, Eq. 34], and the ASM case of (79), in a form 
also involving the ASM bulk statistics u and /x, was obtained by Colomo and Pronko 
Eq. (5.29)], [ , Eq. (3.31)]. 

An expression for Z^^^(l, 1, zi, ^2) = Z^^^ [1,1, zi, Z2) = J2ijlo^n,i,jZ{z^2 i'^ terms of 
a particular Schur function can be obtained by considering (41) with its parameter q set 
to e^*'^/^, and using a formula of Okada [19, Thm. 2.4(1), 2nd equation] for the partition 
function of the six-vertex model with DWBC at g = e^*"^/^. (See Di Francesco and Zinn- 
Justin [17, Eqs. (2.2) & (2.4)] or Biane, Cantini and Sportiello [6, Eq. (1.7)] for further 
details.) It was shown by Biane, Cantini and Sportiello [6, Thm. 1], using this expres- 
sion together with a determinantal identity for certain Schur functions, which is itself ob- 
tained using a general identity involving minors of a matrix, that the numbers An % j satisfy 

deto<iJ<n-l{An,ij) = (-l)"("+l)/2+l {An-lT-'. 

The doubly-refined enumeration of ASMs without bulk statistics, and without reference to 
DPPs, has also been considered in the context of totally symmetric self-complementary plane 
partitions (TSSCPPs). In particular, it was conjectured by Mills, Robbins and Rumsey [46, 
Conj. 3] that the joint distributions of certain pairs of statistics on TSSCPPs in a 2nx2nx2n 
box are equal to the joint distribution of pi and p2 (or ps and p^) on ASM(n) (see also 
Robbins [\ , p. 16]), and this has been proved by Fonseca and Zinn- Justin [26]. Pfaffian and 
constant-term expressions for the associated doubly-refined TSSCPP generating function 
have been obtained by Ishikawa [30, Thm. 1.4, Cor. 7.3 & Cor. 8.2], and, using the result 
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of [26] and Theorem 1, these expressions, as well as certain integral expressions of [26, 
Eqs. (4.9) & (4.14)], also apply to Zf^{l, 1, z,, z^) = ^^^^(1, 1, zi, z^). 

Furthermore, the numbers An,i,j have been considered in the context of a doubly-refined 
Razumov-Stroganov conjecture by Di Francesco [Jf), Sec. 4]. 

Proceeding now to other types of doubly-refined ASM enumeration, for i ^ j and 1 < 
i,i < 4, define Z^^ {x,y, zi, Z2) = J2AeASM{n)^^^^''y^^^^ zl^''^^ Z2^'^\ where it can be checked 
that this is symmetric in i and j. Thus, for the already- considered case, which involves oppo- 
site boundaries of each ASM, Z^^^{x, y, zi, Z2) = Z^'^{x, y, Zi, Z2) = Z^^{x, y, Zi, Z2). It can 
also be checked, using symmetry operations on ASM(n) and (70)-(71), that for the remaining 
cases, which involve adjacent boundaries of each ASM, Z^'^(x, y, zi, Z2) = Z^'^^^x, y, zi, Z2) = 

n(n-l)/2 ^ ^1,4/1 V ± ±^ = ^n{n-l)/2 ^ )n-l ^2,3/1 V ± ^ J^pp gtatistic 

is currently known whose enumerative behavior together with the DPP boundary statis- 
tic pi (or P2) matches that of the statistics associated with two adjacent boundaries of an 
ASM. However, a remarkably simple relation between the generating functions for opposite- 
boundary and adjacent-boundary doubly-refined enumeration of ASMs, for the case without 
bulk statistics (i.e., x = y = 1), has been obtained by Stroganov [ , p. 61]. 

Proceeding to quadruply-refined ASM enumeration, by following an approach similar to 
that of Section 5, but using the form (30), instead of (31), of the Desnanot-Jacobi identity, 
an expression is obtained in [4] for X]AeASM(n) x'^^^^y'^^^^ Zi^^^^ z^^'^^^ z''^'^^^ z^^'^^^ in terms of 
ASM adjacent-boundary doubly- refined generating functions. 

Finally, multiply-refined ASM enumeration involving boundary statistics other than just 
those associated with the first and last row and column of an ASM, but not involving bulk 
statistics, has been studied by Fischer [22, 23], Fischer and Romik [21], and Karklinsky and 
Romik [ v)]. For example, in [22, 33] a relation between opposite-boundary doubly-refined 
ASM enumeration and doubly-refined ASM enumeration involving statistics associated with 
the first and second (or last and second-last) rows (or columns) of an ASM is obtained, 
and a simple formula for the latter is derived, while in [23] a relation between two types of 
ASM triply- refined enumeration, one of which involves any three of the statistics pi, . . . , p^, 
is obtained. No DPP statistic is currently known whose enumerative behavior matches that 
of the additional ASM statistics in these cases. 
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